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We present the scaling properties of A, H, Q, and their anti-particles produced at mid-rapidity in 
Au-\- An collisions at RHIC at ^^snn ~ 200 GcV. The yield of multi-strange baryons per participant 
nucleon increases from peripheral to central collisions more rapidly than the A yield, which appears 
to correspond to an increasing strange quark density of matter produced. The value of the strange 
phase space occupancy factor 73, obtained from a thermal model fit to the data, approaches unity 
for the most central collisions. We also show that the nuclear modification factors, Rcp, of A and 
H are consistent with each other and with that of protons in the transverse momentum range 2.0 
< pt < 5.0 GeV/c. This scaling behaviour is consistent with a scenario of liadron formation from 
constituent quark degrees of freedom through quark recombination or coalescence. 
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Lattice Quantum ChromoDynamics calculations pre- 
dict that a new state of matter, the Quark Gluon Plasma 
(QGP), can be formed at zero baryon density in nu- 
clear collisions when the temperature exceeds 160 — 170 
MeV 0. Strange quarks, whose mass is comparable to 
the critical temperature, are expected to be abundantly 
produced by thermal parton interactions in the high tem- 
perature QGP phase. Due to the corresponding increase 
in the strange quark density, hyperon production is ex- 
pected to be enhanced in high energy nuclear collisions, 
the enhancement increasing with the number of strange 
valence quarks in the hyperon Such an effect has al- 
ready been observed in various fixed-target experiments 
at lower energy by comparing the number of hyperons 
produced per participating nucleon in nucleus-nucleus 
and proton- nucleus collisions In this letter, we 

study the centrality dependence of hyperon production 
in Au + Au collisions at a collision energy of ^/snn — 
200 GeV, which is an order of magnitude higher than 
that previously achieved. We also study the transverse 
momentum dependence of hyperon production in central 
and peripheral collisions in an attempt to shed light upon 
the possible production mechanisms. 

Previous studies have shown that ratios of hadron 
yields in high energy nucleus- nucleus collisions are gen- 
erally well described by statistical models in the grand 
canonical limit, where barvp-chemical potential and tem- 
perature are parameters [a, • A strangeness phase- 
space occupancy factor, 7^, is sometimes introduced to 
describe the extent to which strangeness reaches its equi- 
librium abundance. In this framework, the amount of 
strangeness produced per participating nucleon (Npart) 
is directly related to the value of 7^. The centrality de- 
pendence of 7s therefore provides a quantitative measure 
of strangeness equilibration as a function of system size 
in nucleus-nucleus collisions 9], under the assumption 
that the grand canonical approximation remains valid in 
non-central collisions. 

By contrast, at high transverse momentum, hadrons 
are thought to be produced via incoherent hard scatter- 
ings, which, in the absence of any nuclear medium effects, 
should scale with the number of binary nuclcon-nuclcon 
collisions (Nbinary) 0, 0| • Measurements of hadron 
production in Au + Au collisions at RHIC have shown 
that not only is there a deviation from binary scaling in 
central colhsions 0, 0, but also a distinct difference 
in the scaling behaviour of baryons and mesons in the 
intermediate transverse momentum range, 2 < px < 5 
GeV/c [HQ. A strong particle- type dependence is not 
predicted by conventional Monte Carlo (MC) event sim- 
ulators such as HIJING, where hadron formation in this 
region is dominated by independent parton fragmenta- 
tion jl5l |. On the other hand, quark recombination (co- 



alescence) models have been successful in explaining the 
observed deviation from binary scaling for baryons and 
mesons in central collisions |lall3illaill2|i as well as pro- 
viding an explanation for the particle-type dependence 
of measured azimuthal anisotropics at intermediate pt 
in non-central collisions 'l4|. By extending these previ- 
ous studies to include multi-strange baryons we provide 
a more stringent test of recombination models. Further- 
more, it may allow us to probe the differences between 
strange and light (up and down) quark distributions pro- 
duced in nucleus- nucleus collisions. 

The STAR Time Projection Chamber (TPC) mea- 
sures the trajectories and momenta of charged particles 
produced in each collision in the pseudo-rapidity range 
I77I < 1.8 [23. The detector operates within a solenoidal 
magnetic field of 0.5 Tesla whose axis is aligned with 
the beam. A central trigger barrel, covering the pseudo- 
rapidity region |?7| < 1, and two zero-degree calorime- 
ters are used as trigger detectors. A total of 1.6 x 10^ 
minimum-bias trigger collisions and 1.5 x 10^ central trig- 
ger collisions were used for this analysis. A detailed 
description of the analysis including particle reconstruc- 
tion, track quality, decay vertex topology cuts and cal- 
culation of the detection efficiency can be found else- 
where m m Ei- In this study A(A), S-(S^) and 
(ri ) have been measured in rapidity intervals of 
\y\ < 1, 0.75 and 0.75, respectively. In order to increase 
statistics, the results for ft~ and have been com- 
bined. Within the chosen rapidity intervals the particle 
reconstruction efficiency is a function of transverse mo- 
mentum and lifetime. The efficiency calculations were 
based on the probability of finding Monte Carlo gen- 
erated particles after processing them through a TPC 
detector response simulation, embedding them into real 
events and then reconstructing them as real data. The 
collision centrality was defined by the charged particle 
multiplicity measured in the TPC in the pseudo-rapidity 
range |?7| < 0.5. Five centrality bins were selected cor- 
responding to the following ranges in the total hadronic 
cross section (0 - 5%, 10 - 20%, 20 - 40%, 40 - 60%, 
60 — 80%). The — 5% bin represents the most cen- 
tral collisions and was obtained from the central trig- 
ger sample. The remaining bins were obtained from the 
minimum-bias sample. Due to relatively poor statistics, 
the 5-10% bin and the fl 10-20% and 60-80% bins were 
omitted from this analysis. 

Figure^shows the transverse momentum distributions 
of A(A), S~(S''^) and f2~ + measured at mid-rapidity 
and as function of centrality. The errors shown on the 
data points are statistical only. The A spectra were cor- 
rected for feed-down from multi-strange baryon weak de- 
cays, based upon the measured S and f2 spectra. The 
feed-down correction depends sensitively on both exper- 
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FIG. 1: Transverse momentum distributions of (a) A(A) for 
\y\ < 1.0, (b) H" (H^) for < 0.75 and (c) fi" +11^ for 
y| < 0.75 in Au + Au collisions at ^/snn = 200 GeV as 
a function of centrality. The A spectra were corrected for 
weak decay of H, H° and Q,. Scale factors were applied to 
the spectra for clarity. Only statistical errors are shown. The 
dashed curves show a Boltzmann fit to the A, "ET and fi^ 
data, the fits to the A and are omitted for clarity. 

imental acceptance and the cuts used in the analysis. 
The contribution to the A spectrum from S and 17 de- 
cays is at the 15% level. The feed-down contribution to 
the S spectrum from 17 decays is negligible. The mea- 
sured pt coverage is about 70% for A and 60% for S and 
Q,. The total integrated yields {dN/dy) were extracted 
from Boltzmann fits to the spectra and are presented in 
Table |I| 

The systematic error on the reconstructed yields was 
studied as a function of pt- Three main factors con- 
tribute to the systematic error: (i) subtle differences be- 
tween the Monte Carlo simulation and real data, which 
make the reconstructed yields sensitive to the choice of 
geometric cuts used to improve the signal to background 
ratio, (ii) sensitivity to the method used to subtract the 
remaining background after geometric cuts have been ap- 
plied, and (iii) measured differences in the yield depen- 
dent on the direction of the applied magnetic field. At 
low Pt, the dominant contribution to the systematic er- 
ror is due to the choice of cuts. Here, the systematic 
error was estimated by varying the cuts about their opti- 
mal values and observing the change in the reconstructed 
yield. At high pr the systematic error is dominated by 
the differences observed in the reconstructed yield for the 
two magnetic field settings. In order to determine the 
systematic uncertainty on the total yield and the inverse 
slope parameter for each particle and centrality class, 
Pt dependent systematic errors were added to the data 
points shown in figurenand included in a second fit. The 
systematic errors shown in Table reflect the difference 
between the two fits. We also investigated the choice of 
function used to fit the data. Although the Boltzmann 
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FIG. 2; (a) The corrected integrated yield dN/dy at mid- 
rapidity for H^, A and p divided by Npart, normalized to the 
most peripheral centrality interval (60 — 80%), plotted as a 
function of Npart- The gray, black and dashed bands represent 
the errors on the normalization to the most peripheral bin for 
the H^, A and p. Other errors shown are statistical only, 
(b) 7s as a function of Npart calculated from thermal model 
fits to the measured particle yields (7r,K,p '23], A, H, and 
their anti-particles) at 200 GeV. Values for -|- e~ and p-l-p 
collisions at ^/sNN ~ 91 and 200 GeV respectively and for 
Pb-|-Pb SPS collisions at ^sjvjv ~ 17.2 GeV are shown for 
comparison |2 l l27 l l2^. 

function gave a better fit, an exponential function could 
not be excluded. Exponential fits to the data gave a 5-6% 
higher yield on average and a larger inverse slope param- 
eter by 40-50 MeV. These differences are not included in 
the errors shown in Tabled 

Figure El^a) presents the strange anti-particle yields, 
dN/dy, divided by Npart- For clarity, only statistical er- 
rors are shown. All data points are normalized to the val- 
ues obtained in the most peripheral collisions (centrality 
bin 60 — 80%). The centrality dependence of the anti- 
proton yield is also shown for comparison j24| . Strange 
anti-particles are chosen because all valence quarks must 
have been created in the collision, although similar re- 
sults are also obtained for strange particles. In a geo- 
metrical description of nuclear collisions the number of 
participant nucleons is proportional to the initial overlap- 
ping volume of the colliding nuclei. The integrated yield 
is dominated by the low pt region where particle pro- 
duction originates from mainly soft (non-perturbative) 
processes. The integrated yield per participating nu- 
cleon may be a measure of the formation probability of 
a hadron from the bulk. As such, we would expect it 
to be sensitive to the density of the hadron's constituent 
quarks in the system. We note that there appears to be a 
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Centrality 


0-5% 


10-20% 


20-40% 


40-60% 


60-80% 


< Npart > 


352±3 


235±9 


141±8 


62±9 


21±6 


A 


16.7±0.2±1.1 


10.0±0.1±0.7 


5.53±0.05±0.39 


2.07±0.03±0.14 


0.58±0.01±0.04 




309±1±8 


308±1±8 


303±1±8 


297±2±10 


287±3±10 


A 


12.7±0.2±0.9 


7.7±0.1±0.5 


4.30±0.04±0.30 


1.64±0.03±0.11 


0.48±0.01±0.03 




310±1±7 


309±1±8 


306±1±9 


298±2±10 


282±3±10 




2.17±0.06±0.19 


1.41±0.04±0.08 


0.72±0.02±0.02 


0.26±0.01±0.02 


0.063±0.004±0.003 




335±4±7 


331±4±8 


326±3±6 


325±4±7 


320±8±13 




1.83±0.05±0.20 


1.14±0.04±0.08 


0.62±0.02±0.03 


0.23±0.01±0.02 


0.061±0.004±0.002 




335±4±9 


334±4±9 


327±3±6 


327±5±7 


302±8±16 




0.53±0.04±0.04 




0.17±0.02±0.01 


0.063±0.008±0.004 






353±9±10 




348±15±12 


336±17±13 





TABLE I: Integrated yields dN/dy and inverse slope parameters T (MeV) extracted from a Boltzmann fit to the pr spectra of 
A(A), H~(H)+ and Q,^ + fl^ at mid-rapidity. Statistical and systematic errors are presented. Also shown for each centrality is 
< Npart >, the number of participants, extracted from a Monte Carlo Glauber model calculation ITl|. 



hierachy of particle production dependent upon strange- 
quark content, which has also been observed at lower en- 
ergies 0, Q . This may reflect an increase in the strange 
quark density in more central collisions. 

Thermal-statistical models have been very successful 
in describing particle yields in various systems at differ- 
ent energies Within such models, the densities of 
strange particles, including strange resonances, are gov- 
erned by statistical laws. The possible non-equilibrium 
of strange quarks is taken into account by introducing a 
phase space occupancy factor, 7s. With the measured 
yields of strange baryons and other hadrons, such as 
pions, kaons, protons and their anti-particles we 
have performed a fit using the statistical model described 
in |25j to determine 7^ as a function of the number of 
participants, as shown in figure [^Jb). We find that the 
value of 7s increases from about 0.8 in peripheral col- 
lisions to about 1.0 in central collisions. In each case 
we obtained a freeze-out temperature around 165 MeV. 
According to the model, the A yield depends linearly 
on 7s while the yield of S depends on 7^. This is con- 
sistent with behavior observed in figure HJa). The fact 
that 7s approaches unity when Npart > 150 suggests that 
the strange quark abundance tends to equilibrate as the 
system-size increases. A recent analysis of hadron yields 
in nucleus-nucleus collisions at y^SNN = 17.2 GeV, using 
a different thermal model, also found that 7s approaches 
unity at mid-rapidity in central collisions [26j, whereas 
statistical analyses of elementary e^-l-e" and p+p colli- 
sions at various energies yield a value of 7s significantly 
less than unity [27I Eg . 

We studied the effect of including different combina- 
tions of particles in the fit and found that particle ra- 
tios involving protons and A are important in constrain- 
ing the freeze-out temperature and 7s, respectively. The 
value and centrality dependence of 7s is relatively insen- 



sitive to the inclusion of other particle ratios in the fit. 
The errors shown in figure[5fb) reflect the variation of 7s 
found in this study. 

In order to investigate the scaling behaviour of hy- 
peron production in the intermediate transverse momen- 
tum region, figure|31shows the nuclear modification factor 
(Rcp) ilH for S" -I- S'*' and fl+Tl. The nuclear modi- 
fication factor was found by forming the ratio of the pT 
spectra of the — 5% and 40 — 60% centrality bins, after 
normalising each spectrum to the average number of bi- 
nary collisions, appropriate for each centrality range, ob- 
tained from a Monte Carlo Glauber calculation |l(i llll | . 
The 40 — 60% centrality bin was chosen as the refer- 
ence because of the limited statistics of -|- i7 in the 



Au + ^^^Au at Vs„., = 200 GeV 




h A+ A ^^^1 binary scaling 

0.1 L ^ I ^ I ^ I 

2 4 6 

Transverse momentum p-p (GeV/c) 



FIG. 3: Rcp for -I- H and + ^1 at mid-rapidity 
(centrality interval : — 5% vs. 40 — 60%). A dashed line 
for charged hadrons and gray band for A -I- A are shown as 
comparison. The gray rectangles represent participant and 
binary scahngs. 
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60 — 80% bin. Also shown in figure |2| are the previ- 
ously published results for charged hadrons and A + A 
for the same centrality bins [l^. The dark gray rect- 
angular boxes on the plot represent the expected Rcp 
range for Npart and N^in scalings, indicating the range 
of uncertainty in calculating the number of participants 
and of binary collisions for each centrality. Although 
the pt integrated yield per participating nucleon of S 
increases faster with Npart than for A hyperons, in the 



Republic, FOM and UU of the Netherlands, DAE, DST, 
and CSIR of the Government of India; the Swiss NSF. 



< Pt < 3.5 GeV/c, the pr dependence of 
-I- and fl^ ' 



interval 1 

Rcp for S~ -I- S ' and fi^ -I- ' are similar and coincide 
with the trend previously shown for A -I- A. The Rcp 
of hyperons exhibits little suppression while mesons (ap- 
proximated by the dashed line) have a distinctly different 
trend. The difference in Rcp for baryons and mesons in 
the intermediate px region has previously been discussed 
in the framework of recombination (or coalescence) mod- 
els 0, Q, 0, . The results presented here appear to 
confirm that the difference is dependent upon the number 
of constituent quarks and is not a mass effect. Further 
weight is given to this argument by a recent measure- 
ment of the nuclear modification factor of protons [isj . 
K{892)* [13 and (j) mesons HH. The similarity between 
A and 2 Rcp at intermediate px reinforces the notion of 
a baryon-meson difference. Futhermore, it suggests that 
the strange quark distribution scales with centrality in a 
similar way to up and down quarks, since baryons with 
different strangeness content seem to follow the same pat- 
tern. This observation is consistent with recent elliptic 
flow measurements of A, S and at intermediate |32l |. 

In this letter, we have presented the scaling properties 
of strange baryon production in Au + Au collisions at 
^Jsmn — 200 GeV. By studying the hyperon yields scaled 
by Npart and the centrality dependence of 7s within 
the framework of a thermal model, we have found that 
strangeness equilibrium appears to have been achieved 
in central collisions at RHIC. We have also investigated 
the centrality dependence of the transverse momentum 
distributions of hyperons. We find that hyperon yields 
in central collisions fall below the expectation for binary 
scaling for px > i GeV/c and that the nuclear modifica- 
tion factor Rcp is similar for all hyperons independent 
of their mass or strangeness content. In addition we note 
that the Rcp of hyperons is similar to that of protons, a 
feature that is consistent with models of hadron forma- 
tion based upon quark recombination. 
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